Possible direct method to determine the radius of a star 
from the spectrum of gravitational wave signals II : Spectra for various cases 
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We compute the spectrum and the waveform of gravitational waves generated by the inspiral of a 
disk or a spherical like dust body into a Kerr black hole. We investigate the effect of the radius R 
of the body on gravitational waves and conclude that the radius is inferred from the gravitational 
wave signal irrespective of (1) the form of the body (a disk or a spherical star) (2) the location 
where the shape of the body is determined, (3) the orbital angular momentum of the body, and 
(4) a black hole rotation. We find that when R is much larger than the characteristic length of 
the quasinormal mode frequency, the spectrum has several peaks and the separation of the troughs 
Ao; is proportional to . Thus, we may directly determine the radius of a star in a coalescing 
binary black hole - star system from the observed spectrum of gravitational waves. For example, 
both trough frequency of neutron stars and white dwarfs are within the detectable frequency range 
of some laser interferometers and resonant type detectors so that this effect can be observed in the 
future. We therefore conclude that the spectrum of gravitational waves may provide us important 
signals in gravitational wave astronomy as in optical astronomy. 

PACS number(s): 04.30.Db, 04.25.Nx 
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I. INTRODUCTION 

Researchers in gravitational physics have been ex- 
pecting for the direct detection of gravitational waves 
(GW) using laser interferometers such as LIGO, VIRGO, 
GEO600, TAMA300, and LISA as well as several projects 
using resonant type detectors 0]. The direct detection 
of GW will provide us not only the verification of Gen- 
eral Relativity but also plenty of new aspects in many 
other fields, such as nonlinear physics related to the 
Einstein equations, nuclear and particle physics from 
the estimation of equation of states of the neutron star 
(NS), cosmology from the evidence of inflationary uni- 
verse and from the measurements of cosmological param- 
eters. Therefore, it will open the frontiers of physics, 
gravitational wave physics/astronomy in the 21st cen- 
tury. 

Among many possible sources of GW, coalescing bina- 
ries composed of black holes (BHs) or NSs are probably 
the most promising. From these gravitational wave sig- 
nals, we can trace the astrophysical event. Comparing 
the inspiral waveform with the theoretical templates us- 
ing the matched filtering technique, we may determine 
masses and spins of the BH and star, respectively 0. 
On the other hand, different type of information from 
the inspiral phase may be obtained from the final merg- 
ing phase such as the radius of a star, the equation of 
state of a high density matter, and the structure of the 
stars. Since many physical elements have to be taken into 
account at once to obtain such information in the merg- 
ing phase, little numerical and analytical study have been 
effectively investigated. Recently, Vallisneri M argued a 



possibility to extract the radius of a star from GW. If we 
measure the frequency of GW from the location where 
the star is tidally disrupted, we may find information 
about the central density of the star, and from that, we 
may obtain the radius of a star if we know its equation 
of state. This method as well as other methods using 
the results of numerical simulations, however, may need 
information about the equation of state to determine the 
radius of the star. For example, Vallisneri |^ used the 
result of an incompressible homogeneous Newtonian fluid 
in a Kerr BH g]. 

The first step to tackle this problem of the merging 
phase might be examined by using a BH perturbation ap- 
proach. In this paper, we shall consider GW from either 
an usual star or a NS inspiraling into a BH, emphasizing 
on the size effect of a disk/star on GW using a BH per- 
turbation approach. Many papers along this approach 
have been published on this topic; GW from a test parti- 
cle plunging into a BH, since Sasaki and Nakamura 
reformulated Teukolsky equation Q to make it possible 
to calculate GW induced by a particle plunging into a 
Kerr BH from infinity. Especially, several papers focus 
on phase cancellation effects on GW using this approach. 
In early studies by Nakamura and Sasaki |Q and Haugan, 
Shapiro, and Wasserman GW from a deformed shell 
falling into a Schwarzschild BH were investigated. Saijo, 
Shinkai, and Maeda ||] extended their model to a Kerr 
BH and found that when the central BH is rotating, a 
slightly oblate dust shell minimizes the collapsing energy 
of GW but with a non-zero finite value, which depends 
on the Kerr parameter. 

A BH perturbation approach is not only a toy model to 
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illustrate GW from a binary system. There are two well- 
known facts comparing the results from the BH pertur- 
bation approach with the full numerical simulation. The 
research of the head on collision of two BHs which was 
pioneered by Smarr JlOt and updated by Anninos, Price, 
PuUin, Seidel, and Suen was the first work. They 
compared their numerical results with that of the BH 
perturbation approach ; GW from a test particle falling 
into a Schwarzschild BH by Davis, Ruffini, Press, and 
Price The radiated energy agrees with each other 

within a factor of two. Especially, the waveform of the 
ringing tail has a very accurate coincidence between these 
two cases. The second research was an axisymmetric ro- 
tating star collapse by Stark and Piran |13|. They also 
compared their result with the previous semi-numerical 
result; ring shaped test particles plunging into a Kerr BH 
by Nakamura, Oohara, and Kojima |1J]. Although these 
two models are quite different, the waveforms of these 
two computations give a fairly good agreement. We can 
interpret the reason why these different two results have 
a coincidence as follows. The main body of the rotating 
star first collapses to form a Kerr BH. Then, the remain- 
ing part which is essentially regarded as ring-shaped par- 
ticles fall into the newly formed BH. Therefore, the above 
agreement is not so strange from this physical picture of 
a rotating star collapse. 

These two examples encourage us to apply the BH per- 
turbation approach to the final phase of coalescing binary 
systems. The model worked in this approach may throw 
a milestone on what kind of situation provides us a very 
impressive viewpoint for GW and needs to be followed 
by the detailed, numerical simulations. The result from 
this approach also helps those who work in this field to 
concentrate on such predicted, exciting results. 

Our purpose in this paper is to focus on the size ef- 
fect of a star on GW from a coalescing binary BH - star 
system. We use the BH perturbation approach for the 
sake of the detailed, perspective inve stig ation to find the 
property appears in the GW. In Ref . , we have already 
presented the result with a limited disk with one 

set of parameters. In this paper, we examine the depen- 
dence of the results on (1) the form of the body (a disk 
or a spherical star), (2) the location where the shape of 
the body is determined, (3) the angular momentum of 
the body, and (4) a BH rotation in detail. 

This paper is organized as follows. In Sec. ^ we in- 
troduce the method how to construct the radial wave 
function of GW for disk shaped test particles as well 
as spherical shaped ones using the BH perturbation ap- 
proach. We present our numerical results of GW from 
a disk and a spherical dust star in Sec. III. Section IV 
is devoted to discussion. Throughout this paper, we use 
the geometrized units of G = c = 1 and the metric sign 
as 



II. RADIAL WAVE FUNCTION 

In this section, we explain the method to construct 
the radial wave function of a dust disk and a spherical 
star in Kerr spacetime using the BH perturbation ap- 
proach. Before explaining the construction method of 
the radial wave function, let us briefiy introduce the BH 
perturbation approach Pj5|,p^. This approach is appli- 
cable only if AI ^ /i, where M is the BH mass and 
fi is the particle mass. However, the computational re- 
sult happens to be well-behaved in some cases even if we 
extrapolate the result to /i M as we have already men- 
tioned before. In order to extract GW at infinity with 
sufficient computational accuracy, we have to solve the 
generalized Regge- Wheeler equation (Sasaki-Nakamura 
equation) [Eq. (2.28) of Ref. |]] such as 
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(2.1) 



where r* , Simui{T*), F{r*) and U{r*) are the tortoise co- 
ordinate of a Kerr BH, the source term from a test par- 
ticle of mass fj, [Eq. (2.29) of Ref. B] and two potential 
functions [Eqs. (2.12a) and (2.12b) of Ref. |]], respec- 
tively. To solve Ximui {r* ) using a Green function method, 
we need two independent homogeneous solutions whose 
boundary conditions are given by 
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where k = uj — ma/[2(M + \JM'^ — a^)]. The inhomoge- 
neous solution to Eq. (2.1) becomes 
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where W is the Wronskian, 



(2.2) 



The asymptotic behavior of the radial wave function 
Xi,nLu{r*) is given by 

Xi„,^ir*) = Aimu.e''^'^', (2.3) 

Slmui{r*) in(0) 



Using the amplitude of the radial wave function, the en- 
ergy spectrum and the waveform of GW [Eqs. (3.6) and 
(3.10) of Ref. §] are given by 
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(-00 < iu < 00), (2.4) 
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where r is the coordinate radius from the center of the 
BH, -28^ is spin —2 weighted spherical harmonics, cq 
is a constant given in Ref. 

First, we consider a dust disk plunging into a Kerr BH 
from infinity. We set three assumptions to construct a 
disk in Kerr spacetime. (1) The disk is made up of test 
particles and each of them moves in the equatorial plane 
plunging into the Kerr BH from infinity. (2) The shape 
of the disk is set at the location r = tq. (3) All com- 
ponent particles of the disk have same energy, angular 
momentum, and Carter constant. In this paper, we set 
the disk at the location rg as lOM or 20M, because our 
purpose in this paper is to focus on the merging phase 
of the coalescing binary BH - star . The motion of 
a test particle in Kerr spacetime is written in general as 
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where E — ^lE is the energy, = is the orbital 
angular momentum, C is Cater constant of the particle, 

2Mr + a2. The symbols P, 



S = + cos^ 9, A = r'^ - 
R and 8 are defined by 

P = E{r^ + a^) - aL, 
R^P^- A[r^ + (L, 
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When a particle moves in the equatorial plane or in the 
constant zenith angle. Carter constant is required to be 
zero from the stability of its trajectory. 

Here we explain our effective method to construct the 
radial wave function, instead of computing all the tra- 
jectories of each component particle. We note that the 
geodesic in the equatorial plane of a Kerr BH is charac- 
terized by two parameters; the specific energy [E) and 
the specific angular momentum {Lz). Let t = T{r) and 
(f) — ^{r) express the orbit of the geodesic for given E and 
Lz ■ Since there is time symmetry and azimuthal symme- 
try in Kerr spacetime, we can find another geodesic with 



same energy and angular momentum. Then, we can set 
another geodesic for same E and with shifting time 
and azimuthal angle a.s t — T{r) -t- c* and = $(r) -I- 
where Ct and are constants. After we set the point 
where the representative particle pass through (r = r.i 
and <j) — (pi at t = T{ro)), the orbit of each component 
particle can be expressed as i = T{r) + T{rQ) — T[ri) and 
(f) — c[)(r) + (pi — $(ri). Therefore, it becomes possible to 
set a number of particles with same E and to form 
the disk of radius R whose center is r = rg and </) = $(ro) 
at t = T{rf)). Using this character, we can construct the 
radial wave function only with a function of r . We 
can also choose arbitrary E and Lz for the disk except for 
the condition that each test particle plunges into a Kerr 
BH. The parameter range of Lz in the case of = 1 
which we use in this paper is 

- 2M - \/4Af2 + AMa < Lz < 2M + y/iM^ - 4Ma. 

(2.13) 

We construct the radial wave function of GW for a 
dust disk. Since Kerr spacetime has a property of sym- 
metry for t and 0, we can use this symmetry to construct 
the radial wave function in order to save computational 
time as we explained above. Therefore, the radial wave 
function of the disk can be constructed from that 

of a test particle ^1,'^^*"^'°'' as 
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where /i is the mass of the disk, ro is the location of the 
central point of the disk where we set the circular disk, 
R is the coordinate radius of the disk at r = ro, S* is the 
normalization factor defined by 



S~ I ds — 2 drr<j)o{r), 
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(2.17) 
(2.18) 



In order to construct A\'^^''^' [Eq. ( ^.14 )], we only have 



to prepare t{r), and (/)(r). 

Next, we present how to set a spherical dust star at 
r = tq (ro = lOM, 20M) in order to set a realistic model 
of BH - star merging. We also set three assumptions to 
construct a spherical dust star in Kerr spacetime. (1) 
The star is made up of test particles that move along the 
constant zenith angle of a Kerr BH. (2) The shape of the 
star is set at the location of r = ro, and the center of the 
star is located in the equatorial plane of a Kerr BH. (3) 
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All test particles have the same specific energy, specific 
orbital angular momentum, and Carter constant. From 
the assumptions (1) and (3), we cannot take the orbital 
angular momentum of the particles into consideration. 
Since the orbital angular momentum has almost lost by 
the gravitational wave emission in the previous inspiral 
phase, the zero orbital angular momentum may not be an 
absurd approximation. Our main purpose in this mod- 
eling to construct a spherical star is to confirm that the 
difference of the body (a disk or a spherical star) little 
affect our main conclusion. We expect this proposition is 
true because the size effect in our model does not depend 
on the orbit of the star. We also note that the stability of 
the constant zenith angle trajectory also requires C = 
when we choose E — 1. 

We use the same technique as the dust disk case to 
construct the radial wave function of the spherical dust 
star. We can only adopt this technique to the particles 
in each constant zenith angle, the radial wave function 

becomes a function of 



for the spherical dust star A 
6 as 
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where /i is a mass of the star, tq is the location of the 
central point of the spherical dust star where we set tq — 
lOM, 20M in this paper, i? is a coordinate radius of the 
star, V is the normalization constant of the star defined 
by 



lOM, 20M) than before We place the spectra into 
three categories, classified by the size of the disk. 

We briefiy explain the feature of each type of the spec- 
trum. Type 1 spectrum has the same character as that 
of a test particle; the energy spectrum has one dominant 
characteristic frequency for each m mode where most of 
GW are radiated. This frequency corresponds to the 
quasinormal mode (QNM) of the background BH. The 
real part corresponds to the resonant oscillation of the 
BH, and the imaginary part corresponds to the "damp- 
ing" time scale of GW at the late time. The waveform, 
which has a damping oscillation at the late time, can 
also be described by the QNM frequency at that part. 
Since GW from a disk has a very similar character to 
that from a test particle, we may ignore the size effect 
on GW as far as the diameter of the disk is smaller than 
the characteristic length of the QNM frequency. 

We classify the spectrum to type 2 when the diameter 
of the disk is comparable to the characteristic length of 
the QNM frequency. In this case, the sharp peak near 
the QNM frequency has been weakened. GW still have a 
damping oscillation in the waveform, though the ampli- 
tude of the ringing tail is much smaller than that for a 
test particle due to the phase cancellation effect on GW. 

When the diameter of the disk is larger than the char- 
acteristic length of the QNM frequency, the spectrum has 
several peaks (type 3). One remarkable feature appears 
in the spectrum that the separation of the troughs does 
not largely depend on the orbit. Using this character, 
we find a direct method to determine the radius of a star 
]l5[ . In the waveform, we can hardly recognize the damp- 
ing oscillation at the late time. Therefore, we cannot find 
any characteristic frequency of GW, corresponds to the 
QNM frequency for a test particle in a plunging orbit. 
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Therefore, we have to prepare A'if^^*'^'^^'^\6) and t{r,9), 



(r, 9) for a set of constant zenith angle trajectories in 



order to construct Af^^> [Eq. ( p9| )] 



III. GRAVITATIONAL WAVES FROM THE DISK 
AND THE SPHERICAL STAR 

We show our numerical results for a dust disk plunging 
into a Kerr BH. We calculate GW for a wide range of 
parameters {LJM = 2, 0, -3, a/M ^ 0, 0.5, 0.9, tq = 



A. Dependence on a BH rotation 

First, we discuss the dependence of the energy spectra 
and the waveforms on the radius of the disk and the BH 
rotation, keeping the other set of parameters as rp = 
lOAf, Lz/M — 2. Since we choose Lz/M — 2, we only 
focus on / = m = 2 nonaxisymmetric mode which is 
the dominant mode in all modes in both spectra and 
waveforms. 

For a test particle case, the peak frequency which 
corresponds to the QNM frequency of the BH is dif- 
ferent among these three spectra [Figs. |^ (a), || (a), 
and 1^ (a)] . For example, the spectrum has a peak 
at frequency Muj = 0.33 (Z = to = 2) in Fig. ^ (a) 
(a — 0), frequency Moj = 0.41 (I = m = 2) in Fig. 
(a) [a/M = 0.5), frequency Muj = 0.63 (/ = to = 2) 
in Fig. I (a) (a/M = 0.9), while QNM frequency for 
a = 0, is Muj = 0.37 - 0.089i (/ = 2), for a/M = 0.5 is 
Mij = 0.46 - 0.085i {I ^ m ^ 2), and for a/M = 0.9 is 
Muj = 0.66 - 0M5i {I = m = 2) [Fig. la of Ref. |l|] 
and Fig. 3 (c) of Ref. [p3|]. These results indicate the 



4 



fact that the vibration of the BH has a dominant effect 
on GW for a test particle in a plunging orbit. When we 
turn to look at waveforms, Figs. ^ (a), ^ (a), and ^ (a) 
have damping oscillations at the late time, which can be 
described by the QNM frequency (Fig. 4 of Rcf. |2|] and 
Fig. 4 of Ref. 11). 

For type 1 case, both energy spectra and waveforms 
for a disk have a very similar behavior to the case for 
a test particle. The small size of the disk (for example 
R = 0.785M) does not largely affect energy spectra and 
waveforms. The characteristic behavior in spectra and 
waveforms is almost the same as the case of a test parti- 
cle. We confirm this feature by comparing Fig. |l] (Figs. 
|, |) (a) and (b) for the spectra, and Fig. | (Figs. |, |) 
(a) and (b) for the waveforms. The difference between 
the spectrum from a test particle and that from a small 
size disk appears only in the form factor as we follow 
the construction of ^J^'^'^'' in En. ( 2.14| ). When we look 
at the form factor [Figs. (a), | (a), and |9| (a)], |/2wP 
almost takes the value unity for almost all range of the 
frequency, which is equivalent to the case of a test parti- 
cle. The behavior of the form factor lead the conclusion 
that there is little difference between R = 0.785M and 
R = from the viewpoint of gravitational wave signals 
in the plunging orbit. 

For type 2 case, both energy spectra and waveforms 
have a different feature from those of a test particle. In 
the energy spectrum, the sharp peak near the QNM fre- 
quency has been weakened [Figs. |^ (c), ^ (c), and^ (c)]. 
Since we have defined the range of the radius for type 2 
that the diameter of the disk is comparable to the charac- 
teristic length of the QNM frequency, the range depends 
on the BH rotation. We choose each range for each BH 
rotation for type 2 case as R/M ^ 2 for a/M < 0.5, 
while R/M ~ 1.5 for a/M = 0.9. The difference of the 
radius for each BH rotation is also confirmed from the 
form factor [Figs. @ (b), ^ (b), andjo] (b)] since there is a 
sharp trough in the form factor appears near the QNM 
frequency. When we turn to look at the waveform [Figs. 
@ (c), |5| (c), and ^ (c)], the amplitude near the merging 
phase has been weakened due to the phase cancellation 
effect on GW. In spite of this effect, we can still find a 
damping oscillation quite clearly. 

For type 3 case, both energy spectra and waveforms 
have a completely different feature from the case of type 
1 and type 2; several peaks appear in the spectrum [Figs. 
I (d), I (d), and | (d)]. These peaks also appear in the 
form factor [Figs. (c), || (c), and ^ (c)] and take the 
same behavior to the spectrum. We show the relation 
between the separation frequency of the troughs and the 
radius of the disk in Table |. Remarkably, the relation 
RAlo ^ 1 holds for a wide range of parameters. This 
relation make it possible to measure R directly from the 
observed Aw. We will discuss this point in more detail in 
Sec. JIV} When we turn to look at the waveform [Figs. ^ 
(d), ||7d), and ^ (d)], we can hardly find a damping oscil- 
lation at the late time due to the phase cancellation effect 



on GW. We find that the size effect only appears at the 
merging phase comparing the waveforms with each dif- 
ferent radius. For example, we find the difference within 
the range of -100 < t - r*/M < -30 in Fig. | (a = 0)^^ 
within the range of -100 < t - r* /M < -30 in Fig. | 
{a/M = 0.5), within the range of -130 < t-r*/M < -30 
in Fig. ^ [a/M = 0.9). Although we assume that all the 
component test particles come from infinity, the wave- 
forms in the merging phase might not depend on this as- 
sumption so much. Since the dominant wave is described 
by the QNM frequency for a plunging orbit case, irre- 
spective of the particle energy and angular momentum, 
the phase cancellation effect from several QNM frequency 
waves is responsible to the behavior of the waveforms in 
this merging phase. 



B. Dependence on an orbital angular momentum 

Next, we discuss the dependence of energy spectra and 
waveforms on the radius of the disk and the orbital an- 
gular momentum, keeping the other set of parameters 
as ro = lOM, a/M = 0.9. For a test particle case, the 
spectrum definitely depends on the orbital angular mo- 
mentum [Figs. |(a), (a), [n] (a), and p] (a)]. In fact. 



GW are effectively radiated near the QNM frequency and 
the spectrum takes an additional bump in the case of 
L^/M — 2.63 [Fig. ^ (a)]. Since we do not take into ac- 
count of the radiation reaction force effect in our model, 
the amount of radiated energy diverges when the particle 
takes a circular orbit. Therefore, GW are radiated effec- 
tively when the amount of the orbital angular momentum 
approaches to the value of a circular orbit. In the case 
oi Lz/M = 2, m = 2 mode dominates in I ~ 2 spec- 
trum while for L^/M = —3, m = —2 mode dominates in 
the spectrum within the positive frequency range. This 
tendency turns completely contrary within the negative 
frequency range. For example, m — —2 mode dominates 
in the case of L^/M ~ 2, while m = 2 mode dominates 
for Lz/M ~ —3. There is another feature in the spec- 
trum; each ±m mode has a refiection symmetry at the 
zero frequency. This comes from the fact that the sys- 
tem has an equatorial plane symmetry. When we turn 
to look at waveforms (Figs. ||, |l^, and |l^, some un- 
clear damping oscillation appears in Figs.^J (a) and |l^ 
(a). These behavior might come from the superposition 
of several waves, each of them is described by a single 
QNM frequency. 

The classification of the type, depending on the radius 
of the disk, and the behavior of each type are the same 
as the case of Subsec. HI A , we will briefly describe the 
results in each category. Spectra, waveforms, and form 
factor for type 1 case are shown in Figs.^0 (b),JlC (b), 
g (b), [I (b), in Figs. I (b), PI (b)^ g(b), @ (b), 
and in Figs. | (b), |l6| (b), y (b), |l| (b). Since this 
type little affects the size on GW, the disk can almost 
be regarded as a test particle from the viewpoint of a 
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gravitational wave signal. Spectra, waveforms, and form 
factor for type 2 case are in Figs. @_(c), |l^ (c), |ll| (c), 

P(c), in Figs. (c), |ll (c), (c)71ll (c), and in Figs, 
(c), (c), y (c), m (c). The characteristic QNM 
frequency, as appeared in the spectrum for type 1 case, 
has been weakened due to the phase cancellation effect 
on GW, but we can still look at the dumping oscillation 
in the waveform. Spectra, waveforms, and form factor 
for type 3 case are in Figs. | (d), 
(d), in Figs. Eid), |ll (d), ff 
I (d), (d),1g (d), m (d). 



JlO| (d), (d), ^ 
(d), [l5[(d), and in Figs. 
We can hardly find the 
characteristic frequency in both spectra and waveforms. 
We also show the snapshot of the disk form in Figs. |l^ 
and focusing on the modification of the disk. When 
we turn to look at RAu; = C in Table |, the dependence 
of the orbital angular momentum on C is larger than that 
of the BH rotation although it takes 0.87 < C < 1.16 in 
any case. The value C indeed depends on the orbit, as 
we will explain the structure of C in Sec. |^ in detail. 



C. Dependence on the location where we set the 
circular disk 

Let us consider the effect of the radius and the loca- 
tion, where we define the circular disk, on energy spectra 
and waveforms, keeping the other set of parameters as 
LJM = 2, a/M = 0.9. The circular disk at ro = 20M 
change the shape and when it comes to r = lOM the 
aspect ratio between the radial direction and azimuthal 
angle direction becomes 3 : 1 in the case of Fig. |l|. It 
is, therefore, natural to investigate the dependence on 
tq. We show the spectra, the waveform, and the form 
factor in the case of rg = lOM in Fig s. p[ , ^, and |^ and 
in the case of ro = 20M in Figs. |2[ [24 and p4| . Since 
the spectra with the size range of 1.57 i? < 6.18 have 
several peaks in tq = 20Af case, it may be possible to 
determine the radius of the disk. From Table ||, we have 
RAlo = C - 0.9 for tq = 20Af case. Since C - 1 for 
ro = lOM, the value of C slightly depends on ro. We can 
estimate the radius of the disk only using the separation 
of the peaks appeared in the spectrum. We will explain 
this point in Sec. [V. 

For ro = 20M case (Fig. |4|), there is one remarkable 
difference in the form factor from ro = lOM case (Fig. 

the form factor almost has u = reflection symmetry 
for I = m — 2 mode with the existence of the orbital 
angular momentum. This means that the system is al- 
most axisymmetry. For ro — 20M case, the late stage 



of the shape turns out to be quite axisymmetry. In fact, 
at the end of the evolution, say r ~ 2M, the location of 
the particles approaches to axisymmetric state (Fig. [2TI ). 
This causes less gravitational wave emission and weakens 
the amplitude at the late time. 



D. Spherical star 

Finally, we show energy spectra and waveforms of GW 
from spherical dust stars in the case of type 2 {a/M = 
0.9, R = 1.56M) and type 3 (a/M = 0.9, R = 5.88M). 
We cannot express the radial wave functions as the prod- 
uct of a test particle part and of a form factor part in the 
spherical dust star, because ^l^^^*"^'''^ (6'), t{r,9), (/)(r, 0) 
depends on both r and 9. However, the main property of 
the spectrum from the spherical dust star is the same as 
that from the disk. When the radius of a spherical dust 
star is comparable to the characteristic length L of the 
QNM frequency, a phase cancellation effect on GW ap- 
pears in the spectrum [Fig. 



25 



(a)]. Therefore, the spec- 
trum has a different behavior from a test particle, which 
is the same to the disk case as we mentioned in Subsec. 



Ill A . When the radius of the spherical dust star is much 
larger than L, several peaks appear in the spectrum (Fig. 
^ (b)). The behavior of the spectrum for a spherical star 
is similar to that for a disk case. We also show the de- 
pendence on To in Fig. ^ (b), (c). These results have the 
same behavior to the disk case; the spectrum has smaller 
amplitude for ro = 20M than for ro = IQM . Since all of 
the characters appeared in the spectrum are very similar 
to the disk case, it is also possible for a spherical star 
to determine the radius directly. We summarize the esti- 
mation of the radius in Table [II and |^ using the same 
method as the disk case. The method still works quite 
well, but the relative error rate becomes large when we 
compare the result with disk case. The reason for the 
large error may account for 9 dependence of the orbit. 
Since all of the component particles of the disk have the 
same E and for the disk case, the difference of the or- 
bit comes only from t and <f>. While, we have additional 
9 dependence on the orbit for the dust star case, we can- 
not estimate a simple relationship between R and Auj as 
we will discuss in Sec. IV. Therefore, the relative error 



rate for the dust star case becomes larger than that for 
the dust disk case. In spite of the above reason, RAuj is 
still constant within 10% error. We also show the wave- 
form of GW from the spherical dust star in Fig. The 
waveform is also similar and has the same property to 
the disk case. 



IV. DISCUSSION 

Using a BH perturbation approach, we discuss the size 
effect on GW in the spectra from a dust disk and a 
spherical dust star spiraling into a Kerr BH for a wide 
range of parameters (Lzi spiraling case, < a/M < 0.9, 
ro = lOM, 20M). 

First, we find that when the radius of a disk or a spher- 
ical star is larger than the characteristic length of the 
QNM frequency, the phase cancellation effect from the 
waves generated from each different location of the par- 
ticle appears in the spectrum. Therefore, it is meaningful 
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to compare the diameter of the disk or the spherical star 
with the characteristic length of GW, which is described 
by the QNM frequency for the plunging orbit case. We 
define the characteristic length L that the time delay be- 
tween the earliest particle sX r^ — L and the latest one at 
rg+L is equivalent to the time period of the characteristic 
frequency, 



i(ro + L) - <(ro - L) = 



27r 

ti^QNM 



(4.1) 



Since the shape of the disk changes from the circle to 
a long and narrow shape (Figs. |l^, and |2l|) and 
then swallows into a BH, the definition of L is indeed 
represents the typical length scale of GW. For the disk 
case, the time lag between the earliest and the latest 
component particle plunging into a Kerr BH is Ai = 
t(ro + i?) — t{rQ — R). We find that when R is larger than 
L, the energy spectrum of GW has a different feature 
from that for a test particle. This character comes from 
the phase cancellation effect on GW from the earliest and 
the latest particle of the disk plunging into a Kerr BH. 
We summarize our result in Table 0. For the spherical 
star case, it is rather difficult to argue with such a simple 
interpretation because the orbit has another dependence 
9. In spite of the complicated structure of the geodesic 
for a spherical dust star, C is almost constant in Tables 
[H and IV, and we can interpret them in a similar manner 
to the disk case. We should also argue the bound orbit 
case whether this phase cancellation effect works in the 
inspiral phase of a coalescing binary system. Although 
the characteristic frequency is different from the QNM 
frequency, the basic idea (phase cancellation effects on 
GW) would also be the same in that case. Therefore, 
the rule found in this paper can also apply to the binary 
stars in the inspiral phase and it will be confirmed in the 
future. 

Next, we propose a possibility to determine the radius 



(disk) 



of a star. The energy spectrum [dE / duj)^^-^^ 
from a disk is expressed as 



(particle) 



Irauj 



of GW 



(4.2) 



where (di^/dw)!^'^*"^'''^ is the spectrum from a single test 
particle. The spectrum from a test particle has only one 
peak at the frequency wqnm, so that the square of the 
form factor | fmu) P is responsible for this behavior. The 
existence of several peaks in Figs, [t] (c), || (c), p| (c), [l6| 
(c), (c), ^ (c), ^ (b), (c) can be understood by the 
approximate estimation of f„iu> as 



frr 



(4.3) 



where T' = dT(r) /dr, assuming that only e'"* term de- 
pends on t in Eq. (|2.15| ). The frequency where fmuj takes 



{n + l)7r 
T'R' 



n = 0,1,2, 



(4.4) 



whic h agr ees quite well with the numerical results. Equa- 
tion (4.4) suggests that the separation of peaks in the 
spectrum Alo may be in proportion to R~^. In Table | 
and |l|, we show Alo for various R/M and we find the 
following relation 



R 



C— for i? > L. 

Auj 



(4.5) 



The value C does not depend on the radius and take 
a similar value of ^ 1, but depends on the orbits and 
the initial data, where we set the location of the circular 
disk. In fact, C depends on them within ~ 20%. In the 
physical unit, Av = Alo/[2i:) is given by 



Av = 5kHz 



R 



10km 



= 8Hz 



( 



> —Y 

V 7000km / 



(4.6) 



assuming that C — \. Therefore, for NSs and white 
dwarfs, the frequency band is within the detectable fre- 
quency range of some laser interferometers and resonant 
type detectors 

Finally, we argue from the observational point of view 
whether we have a real situation to determine the radius 
of a tidally disrupted star from GW in our model itself. 
In the final phase of a coalescing binary system, the effect 
of the tidal force and the deformation of the star should 
be taken into account. Carter and Luminet [p3| demon- 
strated the deformation of the star using an affine star 
model with Newtonian particle dynamics (bound orbit). 
They separate the deformation stage into five phase, ac- 
cording to the position of the star, and point out that 
when the star moves inside the Roche radius with the 
condition of /3 = (i?Rocho/-Rp) ^ 3, where i?Roche is 
Roche radius and i?p is a periastron radius of the or- 
bit, the tidal force term rapidly grows and the pressure 
and self gravity terms can be neglected. This leads us to 
the conclusion that we can neglect the effect of the pres- 
sure and the self-gravity when /3 ^ 3. Laguna et al. ]2^ ] 
extended the analysis of Carter and Luminet [ p3| using 
smooth particle hydrodynamics in Schwarzschild space- 
time. They concluded that when (3 ^ 10, we can neglect 
the effect of the self-gravity and the pressure. From these 
previous calculations, we may naturally set the model 
that a star is tidally disrupted by a BH in a certain dis- 
tance and then plunge into a BH. In this case, we may 
neglect the pressure force and the self-gravity effect in 
the first approximation so that each of the component 
fluid clement follows the geodesic. 

To apply our results to the realistic case, three condi- 
tions should be satisfied. (1) The star is tidally disrupted 
at Tdisrupt ^ 6M, which is expressed as 



M 



> 



Psta 



disrupt 



3/J, 

47ri?3^ 



(4.7) 
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where pstar is the density of the star. (2) In order that the 
phase cancehation is effective in the spectrum, R should 
be larger than L. This condition is described as 



t{ro + R)- tiro - R) > T 



27r 

t^QNM 



(4.8) 



(3) The star does not contact to the BH. This condition 
can be written as 



R ^ ^disrupt ^ho 



(4.9) 

= 10km and 
7000km to 



We apply the NS of mass /i = 1.4AfQ with R 

the white dwarf of mass /i = O.SM© with R - 

our model. Then, two conditions [Eqs. (4.7) and (4.8)] 
are rewritten as 



47ri?3 /Tdi, 




-3/2 



M- 



1) 



(^disrupt \ 
M ) 

R{Mujqnm) = M2, 



(4.10) 
(4.11) 



where we use the first order approximation of Tay lor ex- 
pansion to expand the left hand side of Eq. (4.S). And 
the other condition [Eq. ( |4.9| )] gives the lower limit of 
the BH mass as 



M > 



R 



[(- 



M 



0] 



(4.12) 



Taking rdisrupt = 6 A / we t abula te these three condi tions 
[Eqs. dOol ) (|4Tl| ) and ( p^ j in Tables and ^ 



From Tables Vl and VII , we indeed have a possibility to 
determine the radius of the NS for M ~ 2Mq and the 
radius of the white dwarf for M ^ IOOOA/q. 

In the real situation, we should also take into account 
of the pressure and self-gravity effect on GW so that it 
is urgent to confirm our proposal by full 3D numerical 
simulations that we can indeed extract the form factor 
from GW and determine the value of C. In any case, it 
is quite possible that the spectrum of GW may give us 
important information in gravitational wave astronomy 
as in optical astronomy. 
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FIG. 1. Energy spectrum of gravitational waves from a 
disk moving on an equatorial plane in Kerr spacetime whose 

radius is set up at r = lOM in the case of a = 0, Lz/M — 2 
[(a) R = (test particle), (b) R/M = 0.785, (c) R/M = 2.33, 
(d) R/M = 5.88]. Wc only show I = 2 mode. Solid, dashed, 
dash-dotted, dotted, dash-three dotted line denotes the case 
of m = —2, —1, 0, 1, 2, respectively. 



FIG. 2. Energy spectrum of gravitational waves from a 
disk moving on an equatorial plane in Kerr spacetime whose 
radius is set up at r = lOM in the case of a/M = 0.5, 
Lz/M = 2 [(a) i? = (test particle), (b) R/M = 0.785, 
(c) R/M = 2.33, (d) R/M = 5.88]. We only show / = 2 
mode. Solid, dashed, dash-dotted, dotted, dash-three dotted 
line denotes the case of m = —2, —1, 0, 1, 2, respectively. 



FIG. 3. Energy spectrum of gravitational waves from a 
disk moving on an equatorial plane in Kerr spacetime whose 
radius is set up at r = lOAI in the case of a/Af — 0.9, 
Lz/M = 2 [(a) R = (test particle), (b) R/M = 0.785, 
(c) R/M = 1.56, (d) R/M = 5.88]. We only show I = 2 
mode. Solid, dashed, dash-dotted, dotted, dash-three dotted 
line denotes the case of m = —2, —1, 0, 1, 2, respectively. 



FIG. 4. Waveform of gravitational waves from a disk mov- 
ing on an equatorial plane in Kerr spacetime whose radius is 
set up at r = lOM in the case of o = 0, Lz/M = 2 [(a) 
i? = (test particle), (b) R/M = 0.785, (c) R/M = 2.33, (d) 
R/M = 5.88]. We only show I = m = 2 mode, setting the 
observer at the infinity of 6 = 7r/2, ^ = 0. 



FIG. 5. Waveform of gravitational waves from a disk mov- 
ing on an equatorial plane in Kerr spacetime whose radius is 
set up at r = lOM in the case of a/M = 0.5, Lz/M = 2 [(a) 
R = (test particle), (b) R/M = 0.785, (c) R/M = 2.33, (d) 
R/M = 5.88]. We only show I = m = 2 mode, setting the 
observer at the infinity of ^ = 7r/2, ^ = 0. 



FIG. 6. Waveform of gravitational waves from a disk mov- 
ing on an equatorial plane in Kerr spacetime whose radius is 
set up at r = lOM in the case of a/M = 0.9, Lz/M = 2 [(a) 
R = Q (test particle), (b) R/M = 0.785, (c) R/M = 1.56, (d) 
R/M = 5.88]. We only show I = m = 2 mode, setting the 
observer at the infinity of ^ = 7r/2, ^ = 0. 



FIG. 7. Form factor of the dust disk moving on the equa- 
torial plane of a Kerr spacetime whose radius is set up at 

r = lOM in the case of a = 0, Lz/M = 2 [(a) R/M = 0.785, 
(b) R/M = 2.33, (c) R/M = 5.88]. We only show the case 
for Z = m = 2. 
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FIG. 8. Form factor of the dust disk moving on the 
equatorial plane of a Kerr spacetime whose radius is set up 

at r = lOM in the case of a/M = 0.5, L^/M = 2 [(a) 
R/M = 0.785, (b) R/M = 2.33, (c) R/M = 5.88]. We only 
show the case for Z = m = 2. 

FIG. 9. Form factor of the dust disk moving on the 
equatorial plane of a Kerr spacetime whose radius is sot up 
at r = lOM in the case of a/M = 0.9, L^/M = 2 [(a) 
R/M = 0.785, (b) R/M = 1.56, (c) R/M = 5.88]. We only 
show the case for Z = m = 2. 



FIG. 10. Energy spectrum of gravitational waves from a 
disk moving on an equatorial plane in Kerr spacetime whose 
radius is set up at r = lOM in the case of a/M = 0.9, 
L^/M = 2.63 [(a) R = Q (test particle), (b) R/M = 0.785, 
(c) R/M = 1.56, (d) R/M = 5.88]. We only show I = 2 
mode. Solid, dashed, dash-dotted, dotted, dash-three dotted 
line denotes the case of m = —2, —1, 0, 1, 2, respectively. 



FIG. 11. Energy spectrum of gravitational waves from a 
disk moving on an equatorial plane in Kerr spacetime whose 
radius is set up at r = lOM in the case of a/M = 0.9, = 
[(a) R = (test particle), (b) R/M = 0.785, (c) R/M = 1.56, 
(d) R/M = 5.88]. We only show / = 2 mode. Sohd, dashed, 
dash-dotted, dotted, dash-three dotted line denotes the case 
of m = —2, —1, 0, 1, 2, respectively. 



FIG. 12. Energy spectrum of gravitational waves from a 
disk moving on an equatorial plane in Kerr spacetime whose 
radius is set up at r = lOM in the case of a/M = 0.9, 
L,/M = -3 [(a) R = (test particle), (b) R/M = 0.785, 
(c) R/M = 1.56, (d) R/M = 5.88]. We only show I = 2 
mode. Solid, dashed, dash-dotted, dotted, dash-three dotted 
line denotes the CEise of m = —2, —1, 0, 1, 2, respectively. 



FIG. 13. Waveform of gravitational waves from a disk 
moving on an equatorial plane in Kerr spacetime whose radius 
is set up at r = lOM in the case of a/M = 0.9, L^/M = 2.63 
[(a) R = (test particle), (b) R/M = 0.785, (c) R/M = 1.56, 
(d) R/M = 5.88]. We only show I = m = 2 mode, setting the 
observer at the infinity oi 9 = n/2, (j) = 0. 



FIG. 14. Waveform of gravitational waves from a disk 
moving on an equatorial plane in Kerr spacetime whose radius 
is set up at r = lOM in the case of a/M = 0.9, 1/2=0 [(a) 
R = (test particle), (b) R/M = 0.785, (c) R/M = 1.56, (d) 
R/M = 5.88]. We only show I = m = 2 mode, setting the 
observer at the infinity of ^ = 7r/2, = 0. 



FIG. 15. Waveform of gravitational waves from a disk 
moving on an equatorial plane in Kerr spacetime whose radius 

is set up at r = lOM in the case of a/M = 0.9, L^/M = -3 
[(a) R = (test particle), (b) R/M = 0.785, (c) R/M = 1.56, 
(d) R/M = 5.88]. Wc only show I = m = 2 mode, setting the 
observer at the infinity oi d = 71/2, = 0. 

FIG. 16. Form factor of the dust disk moving on the 
equatorial plane of a Kerr spacetime whose radius is set up 
at r = lOM in the case of a/M = 0.9, L^/M = 2.63 [(a) 
R/M = 0.785, (b) R/M = 1.56, (c) R/M = 5.88]. We only 
show the case for Z = m = 2. 

FIG. 17. Form factor of the dust disk moving on the 
equatorial plane of a Kerr spacetime whose radius is set up at 

r = lOM in the case of a/M = 0.9, L^=0 [(a) R/M = 0.785, 
(b) R/M = 1.56, (c) R/M = 5.88]. We only show the case 
for Z = m = 2. 

FIG. 18. Form factor of the dust disk moving on the 
equatorial plane of a Kerr spacetime whose radius is set up 
at r = lOM in the case of a/M = 0.9, L^/M = -3 [(a) 
R/M = 0.785, (b) R/M = 1.56, (c) R/M = 5.88]. We only 
show the case for Z = m = 2. 



FIG. 19. Deformation of the shape of the spherical disk 
whose radius is set up at r = lOM in the case of a/M = 0.9, 
L,/M = 2 [(a) R/M = 1.56, (b) R/M = 5.88]. Solid line 
shows the geodesic for the center of gravity of the disk, while 
circle, square, diamond, and triangle show the edge of the 
disk where the location of the center is at r = lOM, 6M, 4M, 
2M, respectively. 

FIG. 20. Deformation of the shape of the spherical disk 
whose radius is set up at r = lOM in the case of a/Ad — 0.9, 
L^/M = 2.63 [(a) R/M = 1.56, (b) R/M = 5.88]. Solid line 
shows the geodesic for the center of gravity of the disk, while 
circle, square, diamond, and triangle show the edge of the 
disk where the location of the center is at r = lOM, 6M, 4M, 
2M, respectively. 

FIG. 21. Deformation of the shape of the spherical disk 
whose radius is set up at r = 20M in the case of a/M = 0.9, 
L^/M = 2 [(a) R/M = 3.13, (b) R/M = 6.17]. Solid line 
shows the geodesic for the center of gravity of the disk, while 
circle, square, diamond, and triangle show the edge of the disk 
where the location of the center is at r = 20M, lOM, 4M, 
2M, respectively. 
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FIG. 22. Energy spectrum of gravitational waves from a 
disk moving on an equatorial plane in Kerr spacetime whose 

radius is set up at r = 20Af in the case of a/M — 0.9, 
L,/M = 2 [(a) R/M = 1.57, (b) R/M = 6.27]. We only show 
1 = 2 mode. Solid, dashed, dash-dotted, dotted, dash-three 
dotted line denotes the case of m = —2, —1, 0, 1, 2, respec- 
tively. 



FIG. 23. Waveform of gravitational waves from a disk 
moving on an equatorial plane in Kerr spacetime whose radius 
is set up at r = 20M in the case of a/M = 0.9, L^/M = 2 
[(a) R/M = 1.57, (b) R/M = 6.27]. Wc only show / = m = 2 
mode, setting the observer at the infinity oi 9 = n/2, = 0. 



FIG. 24. Form factor of the dust disk moving on the 
equatorial plane of a Kerr spacetime whoso radius is set up 
at r = 20M in the case of a/M = 0.9, L^/M = 2 [(a) 
R/M = 1.57, (b) R/M = 6.27]. We only show I = m = 2 
case. 



FIG. 25. Spectrum of gravitational waves from a star in 
Kerr spacetime in the case of a/M = 0.9 [(a) ro = lOM, 

R/M = 1.56, (b) ro = lOM, R/M = 5.88, (c) n, = 20M, 
R/M = 6.27]. Wc only show I = 2 mode. Solid, dashed, 
dash-dotted, dotted, dash-three dotted line denotes the case 
of m = —2, —1, 0, 1, 2, respectively. 



FIG. 26. Waveform of gravitational waves from a star in 
Kerr spacetime in the case of a/M = 0.9 [(a) ro = lOM, 
R/M = 1.56, (b) ro = lOM, R/M = 5.88, (c) ro = 20M, 
R/M = 6.27]. We only show I = m = 2 mode, setting the 
observer at the infinity of ^ = 7r/2, = 0. 



TABLE I. Comparison with the characteristic length 
which is appeared in energy spectrum of gravitational waves 
to the radius of the disk in the case of ro = lOM. We only 
focus on I = m = 2 mode. 



a/M 


L/M 


R/M 


MAcj 


I/MAlo 


RAuj 








3.09 


0.370 


2.70 


1.14 








3.82 


0.300 


3.33 


1.15 








4.54 


0.255 


3.92 


1.16 








5.22 


0.218 


4.59 


1.14 








5.88 


0.193 


5.18 


1.13 





2 


3.09 


0.340 


2.94 


1.05 





2 


3.83 


0.275 


3.64 


1.05 





2 


4.54 


0.235 


4.26 


1.07 





2 


5.22 


0.203 


4.93 


1.06 





2 


5.88 


0.177 


5.65 


1.04 





-3 


3.09 


0.295 


3.39 


0.911 





-3 


3.83 


0.243 


4.12 


0.929 





-3 


4.54 


0.203 


4.93 


0.921 





-3 


5.22 


0.173 


5.78 


0.904 





-3 


5.88 


0.148 


6.76 


0.870 


0.5 





3.09 


0.340 


2.94 


1.05 


0.5 





3.83 


0.275 


3.64 


1.05 


0.5 





4.54 


0.235 


4.26 


1.07 


0.5 





5.22 


0.203 


4.93 


1.06 


0.5 





5.88 


0.177 


5.65 


1.04 


0.5 


2 


3.09 


0.335 


2.99 


1.03 


0.5 


2 


3.83 


0.270 


3.70 


1.03 


0.5 


2 


4.54 


0.233 


4.29 


1.06 


0.5 


2 


5.22 


0.198 


5.05 


1.03 


0.5 


2 


5.88 


0.177 


5.65 


1.04 


0.5 


-3 


3.09 


0.300 


3.33 


0.927 


0.5 


-3 


3.83 


0.248 


4.03 


0.949 


0.5 


-3 


4.54 


0.210 


4.76 


0.953 


0.5 


-3 


5.22 


0.177 


5.65 


0.924 


0.5 


-3 


5.88 


0.160 


6.25 


0.940 


0.9 





3.09 


0.370 


2.70 


1.14 


0.9 





3.83 


0.300 


3.33 


1.15 


0.9 





4.54 


0.255 


3.92 


1.16 


0.9 





5.22 


0.220 


4.55 


1.15 


0.9 





5.88 


0.197 


5.08 


1.16 


0.9 


2 


3.09 


0.335 


2.99 


1.03 


0.9 


2 


3.83 


0.270 


3.70 


1.03 


0.9 


2 


4.54 


0.233 


4.29 


1.06 


0.9 


2 


5.22 


0.200 


5.00 


1.04 


0.9 


2 


5.88 


0.177 


5.65 


1.04 


0.9 


-3 


3.09 


0.310 


3.23 


0.957 


0.9 


-3 


3.83 


0.250 


4.00 


0.956 


0.9 


-3 


4.54 


0.213 


4.69 


0.967 


0.9 


-3 


5.22 


0.183 


5.46 


0.956 


0.9 


-3 


5.88 


0.163 


6.13 


0.958 
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TABLE II. Comparison with the characteristic length 
which is appeared in energy spectrum of gravitational waves 
to the radius of the disk in the case of ro = 20M. We only 
focus on Z = m = 2 mode. 



TABLE IV. Comparison with the characteristic length 
which is appeared in energy spectrum of gravitational waves 
to the radius of the star in the case of ro — 20M. We only 
focus on Z = m = 2 mode. 



a/M 


L/M 


R/M 


MAuj 


l/MAu 


RAlj 


a/M 


R/M 


MAuj 


1/MAlo RAlo 








3.13 


0.300 


3.33 


0.938 





6.18 


0.260 


3.85 


1.61 








4.67 


0.192 


5.21 


0.896 


0.5 


4.67 


0.285 


3.51 


1.33 








6.18 


0.154 


6.49 


0.951 


0.5 


6.18 


0.250 


4.00 


1.54 





2 


3.13 


0.290 


3.45 


0.907 


0.9 


4.67 


0.260 


3.85 


1.21 





2 


4.67 


0.182 


5.49 


0.849 


0.9 


6.18 


0.245 


4.08 


1.51 





2 


6.18 


0.146 


6.85 


0.902 















—3 


3.13 


0.270 


3.70 


0.844 















—3 


4.67 


0.172 


5.81 


0.803 















—3 


6.18 


0.130 


7.69 


0.803 












0.5 





3.13 


0.300 


3.33 


0.938 


TABLE V. Characteristic length of gravitational waves 


0.5 





4.67 


0.188 


5.32 


0.877 


for a 


plunging orbit. 


We define the characteristic length L 


0.5 





6.18 


0.141 


7.09 


0.871 


with satisfying Eq. (k 


w. We solve geodesic equation numer- 


0.5 


2 


3.14 


0.290 


3.45 


0.907 


ically 


in order to find the characteristic length. We only use 


0.5 


2 


4.67 


0.183 


5.46 


0.854 


I = m 


= 2 mode of the QNM frequency. 




0.5 


2 


6.18 


0.138 


7.25 


0.852 












0.5 


-3 


3.14 


0.275 


3.64 


0.860 


a/M 






iro = 10A/ 


I/ro=20M 


0.5 


-3 


4.71 


0.172 


5.81 


0.803 










2.78 


2.31 


0.5 


-3 


6.27 


0.127 


7.87 


0.785 





2 




3.03 


2.42 


0.9 





3.13 


0.300 


3.33 


0.938 





-3 




2.42 


2.16 


0.9 





4.67 


0.195 


5.13 


0.910 


0.5 







2.22 


1.85 


0.9 





6.17 


0.141 


7.09 


0.871 


0.5 


2 




2.44 


1.94 


0.9 


2 


3.13 


0.285 


3.51 


0.891 


0.5 


-3 




1.99 


1.74 


0.9 


2 


4.67 


0.185 


5.41 


0.863 


0.9 







1.54 


1.29 


0.9 


2 


6.17 


0.133 


7.52 


0.822 


0.9 


2 




1.71 


1.36 


0.9 


-3 


3.13 


0.275 


3.64 


0.860 


0.9 


-3 




1.42 


1.22 


0.9 


-3 


4.67 


0.173 


5.78 


0.807 












0.9 


-3 


6.17 


0.150 


6.67 


0.927 













TABLE III. Comparison with the characteristic length 
which is appeared in energy spectrum of gravitational waves 
to the radius of the star in the case of ro — 10A4. We only 
focus on Z = m = 2 mode. 



TABLE VI. Allowed mass of the BH (Eqs. and 
( [l.ll[ )) for a NS (/i = I.IM©, R = 10km) to determine the 
radius of the star. A4i and M2 de note the upper limit of the 
mass from Eqs. (4.10) and (tl.lll) and M3 denote the lower 
limit of the mass from Eq. (|4.12|), respectively. 



a/M 


R/M 


MAcj 


1/MAlu 


RAuj 





5.22 


0.175 


5.71 


0.914 





5.88 


0.145 


6.90 


0.852 


0.5 


5.22 


0.175 


5.71 


0.914 


0.5 


5.88 


0.145 


6.90 


0.852 


0.9 


5.22 


0.170 


5.88 


0.888 


0.9 


5.88 


0.150 


6.67 


0.881 



a/M 


L,/M 


Mi/Mq 


M2/M0 


M3/M0 








2.07 


2.07 


1.69 





2 


2.07 


2.06 


1.69 





-3 


2.07 


2.05 


1.69 


0.5 





2.07 


2.55 


1.64 


0.5 


2 


2.07 


2.75 


1.64 


0.5 


-3 


2.07 


2.28 


1.64 


0.9 





2.07 


3.56 


1.48 


0.9 


2 


2.07 


4.12 


1.48 


0.9 


-3 


2.07 


3.00 


1.48 
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TABLE VII. Allowed mass of BH (Eqs. ( ^.lOj ) and ( [4.1l[ )) 
for a white dwarf {jj, = Q.5Mq, R = 7,000km) to determine 
the radius of the star. Mi and M2 de note the upper limit 
of the mass from Eqs. (i.lO) and ([4.1 l|) and M3 denote the 

.121), respectively. 



( p~ro| ) and 
lower limit of the mass from Eq. 



a/M 


L,/M 


Ml /Mo 


M2/MQ 


M3/M0 








64300 


1450 


1180 





2 


64300 


1440 


1180 





-3 


64300 


1440 


1180 


0.5 





64300 


1780 


1150 


0.5 


2 


64300 


1920 


1150 


0.5 


-3 


64300 


1600 


1150 


0.9 





64300 


2490 


1040 


0.9 


2 


64300 


2880 


1040 


0.9 


-3 


64300 


2100 


1040 
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